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Abstract 

An exact renormalization group equation is derived for the free energy 
of matrix models. The renormalization group equation turns out to be 
nonlinear for matrix models, as opposed to linear for vector models. An 
algorithm for determining the critical coupling constant and the critical 
exponent is obtained. As concrete examples, one-matrix models with 
one and two coupling constants are analyzed and the exact values of the 
critical coupling constant and the associated critical exponent are found. 
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The matrix model provides the possibihty for a nonperturbative treatment of 
two-dimensional quantum gravity Exact solutions of the matrix model have 

been obtained for two-dimensional quantum gravity coupled to minimal conformal 
matter with central charge c < 1. It is important to solve matrix model representa- 
tions of two-dimensional quantum gravity coupled to conformal matter with central 
charge c > 1. They are interesting not only as statistical systems on a randomly 
triangulated surface, but also as string models in a target space of arbitrary dimen- 
sions. Although one can easily write down matrix model candidates for cases with 
c > 1 0, |], they are so far not solvable. 

As an approximation scheme for obtaining critical coupling constants and critical 
exponents for such unsolved matrix models, Brezin and Zinn- Justin have proposed 
a renormalization group approach |0. A similar approach has been advocated pre- 
viously for the expansion in various other contexts 0. Several groups have 
examined the consequences of such an approach P], suggesting that the result for 
universal quantities (critical exponents) does not improve even if one evaluates beta 
functions perturbatively up to the first few orders. Therefore in order to make full 
use of such a scheme, we first need to confirm that the approximation method pro- 
duces correct results for the exactly solved cases. The 0{N) vector model, describing 
a discretized one-dimensional quantum gravity, is known to share its double-scaling 
properties with the matrix model for a discretized two-dimensional quantum gravity 
0. Recently we have analyzed the vector model by means of the renormalization 
group approach and have clarified its validity and meaning [|T^]. As with the matrix 
models, a perturbative evaluation of critical exponents for the vector model does 
not improve up to a few orders. Instead of a perturbative approach, we have ob- 
tained an exact difference equation relating the free energy for neighbouring values 
of A^, with the coupling constant shifts of order 1/A^ in infinitely many coupling 
constants. Our crucial observation was that there hold infinitely many identities 
which express the freedom to reparametrize the field space. Owing to these identi- 
ties, we can re-express the flow in the infinite dimensional coupling constant space 
as an effective flow in a finite dimensional effective coupling constant space. Critical 
coupling constants and critical exponents arise as fixed points of the effective beta 
function, and eigenvalues of the derivative matrix of the beta function respectively. 
We can understand the reason why the perturbative evaluation of beta functions 
was inadequate at least in the vector model; we have to collect all contributions 
from higher induced couplings by using the reparametrization identities to obtain a 
beta function. 

The purpose of this paper is to derive an exact renormalization group equation 
for matrix models and examine its validity. Using the eigenvalue representation, we 
find a difference equation relating the partition function Zi^^i{g) to the partition 
function Z^lg + 5g) with the coupling constant shifts 5gk of order 1/A^ in infinitely 
many coupling constants. We find that the induced interactions contain terms which 
are absent in usual matrix models, at least for c < 1. As with the vector model, 
we can derive infinitely many identities expressing the freedom to reparametrize the 
matrix variables. These identities enable us to reduce the renormalization group 
equation to an effective renormalization group equation with a finite number of 
coupling constants. We observe that the effective renormalization group equation 
is nonlinear in derivatives of the free energy with respect to coupling constants. In 



spite of this nonlinearity, we can provide an algorithm to compute fixed points and 
critical exponents. To illustrate this technique, we explicitly obtain fixed points and 
critical exponents together with the operator content for the one-matrix model with 
one and two coupling constants. These agree completely with the exact result. 

The partition function Z^ig) of the matrix model is defined by an integral over 
an X hermitian matrix $ with a generic polynomial potential V^($) 

Z^{g) = j f/^^<|.exp [-NivVm , V{^) = f (1) 

k=l 

The cubic interaction with a single coupling constant g corresponds to V{(^) = 
1$^ + 1$^. We can integrate over the angular variables to obtain an integral over 
the eigenvalues {Xj} [|| 



N 



N 



-Nj2V{Xi 



i=l 



(2) 



Z^ig) = Cn Y[dXi A%{{Xj}) exp 
•' 1=1 

where Ajy denotes the Van der Monde determinant A7v({Aj}) = ni<i<j<Ar('^i ~ ^j), 
andc^ = 7r^(^-i)/Vn;Lip! • 

In order to relate Z^+i to Z^, we shall integrate the {N + l)-th eigenvalue A in 

Zn+1 

ZN+iig) = 1^(^+1)^$ e-(^+^)*^^(*) 

= CN+1 / A^({A,}) e-(^+^)S.^i^(^^) /rfAll |A - A.^ e-(^+^)^W 

i=i •' 1=1 

= '^^+1 f (1^2 ^ ^-{N+l)tT f ^-{N+l)V{X)+2trlog\X-'l>\_ ^g^j 

Cn J J 

The A integral can be evaluated by a saddle point method as a power series in 
around the saddle point, since the effective potential (A^ + l)V^(A) — 2tr log |A — $| 
is of order 0{N^). The saddle point equation 



2 1 2 °° tr <J)" 



determines the saddle point A*^ as a linear combination of not only tr but also non- 
linear terms like tr$"tr$™, which are not present in the original model. Therefore 
by inserting the saddle point X^ into the partition function, we find that infinitely 
many operators of the form tr tr . . . are induced 

ZN+l{g) Cn+1 



ZN{g) Cn 



exp 



- tr V{^) - NV{X') + 2 tr log I A'' - $1 + 0(A^") ) (5) 



from the renormalization transformation with respect to A^. Here the normalized 
average with respect to the measure exp [— A?^ tr \^($)] is denoted by ( ). 

The situation, however, simplifies if we appeal to the large- A^ limit. In this 
limit a multi-point function of ?7(A^)-invariant operators O, O', ■ ■ ■ factorizes into a 
product of one-point functions 

{00'---) = {0){0')--- + 0l^-^y (6) 



In view of this factorization property, eq.(^ can be written as 



ZN{g) 



Cn+1 

Cn 



exp 



(tr\/(<l>)) -iV1/(A) + 2(tr log |A - $1) + 0(A^' 



(7) 



where A is a function of (tr^"^) defined by averaging the saddle point equation 
eq.(i 

\N A-$/ ^oA"+i\A^ / 



V'CX) 



(8) 



Let us introduce the free energy which is normalized to vanish for the Gaussian 
model 



F{N,g)^-—\o^ 



(9) 



ZN{gi = 0,g2 = l,gk = {k > 3)) 

where Zj^ig^ = 0, ^2 = 1, ^fc = (A; > 3)) = 2^/^{7i/N)^''/'^. By taking the ^ cx) 
limit, we find the following differential equation as a renormalization group equation 
for the free energy 



V(A)-2(^tr loj 



A-$ 



3 ^ / 1 



(10) 

We observe that this equation describes a "flow" in an infinite dimensional coupling 
constant space which is enlarged to include interactions like tr$™tr$". We em- 
phasize that this renormalization group equation does not involve a perturbation 
with respect to the coupling constants, unlike the approximation schemes proposed 
in ref . . 

We should be careful in identifying the renormalization group flow in the coupling 
constant space, since all the correlators of tr appearing in the right hand side 
of eq.(|TU]) are not independent. In fact, the partition function Z^i^g) should be 
invariant under reparametrization of matrix variables regular at the origin, $ — > 
$ + e^^"^^^ {n > —1). Then we can construct a tower of identities expressing this 
reparametrization freedom of matrix variables, as with the case of the vector model 



10| . The reparametrization identities can be obtained from the usual procedure to 



derive the Schwinger- Dyson equation |Tl 



J d^'$tr^(<l>"+^exp [-iVtry($)]) = 0, 



(n > -11 



:iii 



The reparametrization identities naturally form the Virasoro vacuum conditions on 
the partition function. After using the factorization property (^ in the large-A^ 
limit, they become 



^^i_tr <l>*\ /4tr 



A^ 



^i_tr (<I>"+V(<I>) 



{n > -1). 



(12) 



Since the one-point function is nothing but the derivative of the free energy with 
respect to coupling constants, (tr$'^)/A^ = kdF/dgk, the reparametrization iden- 
tities provide infinitely many relations among derivatives with respect to different 
coupling constants. 



In the case of the vector model, we have shown |TD[ that one can ehminate 
redundant one-point functions in the right hand side of eq.(0) except for the first 
few ones corresponding to the original potential, by applying the reparametrization 
identities recursively. To facilitate the procedure, it is convenient to define the 
expectation value of the resolvent W{z) = (tr(l/(2; — $)))/A^, with respect to 
which the reparametrization identities ( [T^ ) take the form ^ 

W{zy - V\z)Wiz) + Q{z- V) = 0, Q{z; V)^^. ^^^^-^ " ^)'"' 

(13) 

If we are to begin with the m-th order potential, Q{z; V) is a polynomial in ( tr ^^)/N 
for 1 < < m — 2. Let us choose a coupling constant subspace: gi = 0,g2 = 1, gk = 
{k > m + 1) with gs, - ■ ■ ,gm arbitrary. The first two equations of ([T^) 



= (^tr$) + E^?^(^tr l = (^tr$2) + f ..(Itr*'^) (14) 

enable us to express (tr$)/A^ and (tr$^)/A^ in terms of (tr$'^)/A^ for 3 < k < 
m. Therefore the reparametrization identities allow us to express all of the higher 
induced interactions in terms of ^tr^'^^ /N = kdF/dgk for 3 < < m, and in this 
sense the theory is renormalizable in the large- limit. 

To summarize, the complete set of nonlinear renormalization group equation for 
the one-matrix model consists of the following three equations. The reparametriza- 
tion identities become the so-called loop equation 

W{zf - V\Z)W{Z) + Q . . . , . . = 0. (15) 

V og^ dgmj 

The reparametrization identities simplify the saddle point equation into 



V'{X) -4Qix-gs,---,gm;—,---,-^j = 0. 



(16) 



Finally the effective renormalization group equation becomes 



G («. - -1 - E — + V(A) - 2 log A - 2 £^ (W(.) - -) (.17) 



The above equations show that the renormalization group equation for the matrix 
model is inevitably nonlinear with respect to dF{g) / dgu- Nonlinearity emerges 
out of the renormalization group transformation (induced interactions of the form 
tr $™ tr ■ ■ ■) and from the procedure to eliminate redundant one-point functions. 

Now we shall give an algorithm to determine fixed points and critical exponents 
as well as the free energy F = Y^hLo ^'"^^ fhid) ■ For simplicity let us consider the 
case of the single cubic coupling for the one-matrix model. We first restrict ourselves 



to the sphere topology. The right hand side of the renormahzation group equation 
has nonhnear terms {g = Qs) 

G(«.f)^|/..(.)(f)". (18) 

If we expand the coefficients in powers of the couphng constant around a fixed 
point g^, they are regular 

oo 

Pn{g) = T.Pnk{g-g*)'- (19) 

fc=0 

The non-analyticity of the free energy should stem from solving the differential 
equation (|1^. We assume that the free energy consists of singular and regular 
terms 

F(iV,(7) = F-^(iV,^?) + F-"3(iV,(7), 

oo oo 

F-^ = ^iv-2vr(^), fr{g) = j:<^k{g-g.r, (20) 

h=0 k=0 



00 00 



h=0 k=0 

By definition, the coefficient 60 7^ 0. By comparing the power series expansion of 
the renormahzation group equation, we find consistency conditions for the above 
expansion to be valid. The most singular term {g — g*y~^ determines the fixed 
point 

00 

= 6o7E/^no^<"'- (22) 

n=l 

The next singular term {g — g^)^ determines the critical exponent 7 since &o 7^ 



260 = &o7 



I3nina1 ^ + 2a2 5Z Pnan{n - !)<" 



.n=l n=2 



(23) 



They contain two coefficients ai and 02 of the regular part. However, the consistency 
conditions for terms of order {g — g*Y and {g — g*Y determine ai and 02 respectively 

00 cx) 00 00 

2ai = J2 Pnial 2a2 = J2 /^n2< + 2a2 E Pnina1-' + 2a2^ ^ /3„on(n-l)ar2. (24) 

n=0 n=0 n=l n=2 



In terms of G, eqs. (|2^ ) - (|2^) are equivalent to 
2 

= G^ai{g*,ai), - = G^ga^{g^,ai) +2a2 G^aiai{g*,ai) , (25) 
7 

2ai = G^g{g^,ai), 2a2 = ^G^gg{g^,ai) + 2a2 G^ga^{g^,ai) + 2a2^ G^aiaAg*,ai). 

By solving these four equations, we can determine four quantities, namely the fixed 
point the coefficient oi, the critical exponent 7, and the coefficient 02. The 
susceptibility exponent 70 for the sphere topology is related to 7 via 7 = 2 — 70. 



In the one coupling case at hand, G{g, ai) is given by 



G{g, ai) = -1 - |ai + + ^A^^ - 2 log A 



(26) 



/ dz z + gz^ — \/ (z + gz"^)^ — A{1 + gz — g"^ + ^^g^o-i) 

J±cx> I 



with A = X{g,ai) being the zero of the square root. We have found that the set of 
equations (^5]) have three solutions; two of them are the nontrivial ultraviolet fixed 
point which agree with the exact result for the m = 2 critical points 

= ±(432)-^/^ 2-70 = 5/2, ai = ±3-^/^2(1 - ^3), aa = -4(1773 + 30). (27) 

The other fixed point is located at the origin: g* = and 2 — 70 = —4 which agrees 
with the naive scaling as expected. At each fixed point, the other consistency con- 
ditions determine all the other coefficients {k > 0) and bk/bo [k > 1) recursively 
except for the overall normalization of the singular term bo. Although the equation 
for 02 is quadratic, we choose the branch which continues to the unique solution at 
the origin (the Gaussian fixed point). In this way we can obtain the series expan- 
sion of the sphere free energy around the fixed point, up to one integration constant. 
This situation is completely analogous to the vector model case. 

Up to now we have implicitly assumed that the exponent 7 is an irrational num- 



ber. In the case of the rational exponent as in eq. (p7D , the singular part is still well 
defined, but its higher power will contribute to terms of integer powers of {g — g*) 
[]. Therefore the equation to determine the regular part acquires a contribution from 
the singular part at higher orders and hence equations to determine the regular and 
singular parts mix each other. This mixing does not occur in eqs. (p^) which deter- 
mine the coefficients Ci and 02 at least for the solutions (|27|) . Since susceptibility 
exponents for c < 1 are known to be 7 > 2 for exact solutions, eqs. (P^ are also 
valid in these c < 1 cases. 

It is straightforward to generalize the algorithm to the case of many coupling 
constants. The expansion around the fixed point g^ is given in terms of Agi = gi—gi^. 

00 

frig) = E«^"""(^?*)^^?n---A^7,„ 

n=0 

00 

frig) = j:iy^^g^r'T.b''-'-ig*)^gn---Ag.^, (28) 



n=0 



where we have used a transformation matrix Vi-' which diagonalizes the exponent 
matrix. We obtain by expanding the renormalization group equation in powers of 

da^ ' 7(fc) daWa''' 

2a^ = ^ , 2a^^ = i,^ + a^'=-^ + a^'=-^ + 2aV-^, (29) 
ogi 2 ogidgj ogjOa'^ ogiOa^ oa^oa'^ 



^ Since the distinction between the regular and the singular parts of the free energy becomes 
meaningless for integer 7, we assume that 7 is not an integer. 



where they come from terms of order O (^{Vi^ Agj)^(^)^^j , O (^{Vi''Agk)^(^)-\Vj''Ag, 
0{Agi) and 0{AgiAgj), respectively. These four form a closed set of equations to 
determine fixed points and critical exponents as in the one coupling case. 

As an illuminating example, we analyze the case of two couplings {g^ and g^) for 
the one- matrix model. We find five solutions for the set of equations (^9]) . Firstly, we 
find fixed points at ((73*, (74*) = (±0.3066 . . . , 0.02532 . . .). The free energy behaves 
in their neighbourhoods as 

F^'-^{N,gs,g,) ~ b-[Ags±i-AA72 . . .yAg^^/'+b' ■[Ag,±{-6.209 . . .yAg^j'/^ (30) 

The positions of the critical points and the transformation matrices V^/ in eq. (|28|) can 
be expressed as roots of algebraic equations, though we present them in numerical 
forms for brevity in eq. (|30D . We find two singular terms with the critical exponents 
7 = 7/3 and 7/2 at this fixed point. They correspond to the gravitational dressing of 
operators with bare conformal dimension Aq = —1/5 and 0, respectively [0. This 
fixed point should describe the Lee- Yang edge singularity. The fixed point and the 
critical exponents agree with the exact result for the m = 3 critical point. We also 
have the m = 2 fixed points which we have just seen in the one coupling constant 
case eq.(^), ((73*, (74*) = (±(432)^^/"^, 0). The critical behaviors around these points 
are 



F^'-%N,g,,g,) 



Ag3±i- ^1^9^ 



5/2 



+ b'-[Ag,]-'. (31) 



The first critical exponent 5/2 corresponds to the gravitationally dressed cosmo- 
logical term for the m = 2 critical point. We find no other positive exponent, in 
agreement with the fact that there is no primary field other than the identity in 
the matter conformal field theory before gravitational dressing. The last one is the 
trivial fixed point at ((73*, (74*) = (0,0) which gives the naive scaling behavior 

i^""'(iV, ^3, 9d-b- [Ag,]-' + b' ■ [Ag,]-\ (32) 



In the exact solution it is known that there are the m = 3 critical points and 
the m = 2 critical lines in the two dimensional coupling constant space. We observe 
that the critical behavior of the second term in eqs. (pO|)-(|3^) is realized at each 
fixed point when one approaches the fixed point along the m = 2 critical line. It is 
interesting that all of the above solutions that we find turn out to be on the m = 2 
critical lines. However, let us note that not every point on the m = 2 critical lines 
is a fixed point, as with the case of the vector model. 

If we compare our result with the exact one, another solution corresponding to 
the m = 2 critical point ((73*, (74*) = (0, —1/12) is expected to be present. It is likely 
that such a solution exists, though we have not proved it due to the subtlety arising 
from the Z2-invariance of the critical potential. 

Next we turn to the scaling behavior of the higher genus contributions fh{g) 
to the free energy F = Y^'h=oN~'^^ fh{g). Again for brevity we consider the single 
coupling case. We can start from the difference equation and take into account 
the shift 5gk of higher orders in 1/A^. A differential equation with additional terms 
can be derived from the difference equation by retaining the higher order terms in 




expanding in If we expand the free energy into series, we can separate the 
partial differential equation into a set of ordinary differential equations for each genus 
contribution. It is important to realize that the additional contributions introduced 
into the right hand side G of the renormalization group equation carry additional 
powers of Therefore they will not appear in the coefficient I3n{g) 

F{N, g) = rig, N) + N~'^ ^ ^n{g) E 11 (33) 

h=0 n=l hi + --- + hn = h, i=l "'3 

The inhomogeneous term r{g, N) does get various higher order contributions includ- 
ing those terms from fh' for h' < h — 1. Therefore we obtain 

oo n J£ 

2il~h)Mg) = r,{g) + j:f3n{g) E (34) 

hj > 

where rh{g) is defined by r{g,N) = N^'^^rh[g). Only the inhomogeneous term 
depends on genus, whereas the beta functions I3n{g) for n > 1 are universal for any 
genus. Once we realize this structure of the renormalization group equation, we can 
repeat the same argument as in the case of the sphere to obtain the scaling behavior 
of the higher genus contributions. We find the fixed point condition is the same as 
eq. (p2D except that the coefficient 6o is now replaced by of the singular term of 
the genus h free energy. The critical exponent condition becomes 



2(1 - h)h^, = fooS'^ 



E f^nina^' + 2a2 E Mn - l^^^ 

.n=l n=2 



(35) 



We find that the coefficients ai and 02 of the regular term are needed to fix these 
equations. However, it is notable that they are precisely those regular terms for 
the sphere free energy. Therefore the same conditions as the sphere case (p^ ) are 
sufficient to determine them. In this way we find that the fixed point is universal 
for any genus and the critical exponent 7'' for genus h is given by 

7^^ = (1 - /i)7i, 71 + 70 = 2. (36) 

This result explains the double scaling behavior for the singular part of the free 
energy 

00 00 

F-"^(iV,^) = J2N-'^fng) = E N-'\g - g*f-'''-''X 

h=0 h=0 

= {g- g.f-''^n%N^/^^{g - (?.)). (37) 



So far we have employed the eigenvalue representation of the matrix integral and 
have attempted to integrate the (A^ + l)-th eigenvalue. We can, however, also use 
the coset technique proposed in ref.0 where one integrates over the (A^ + l)-th row 
and column vector $Ar+ij = $]^^_i, (1 < j < A^) and the singlet $7v+iAr+i = a 
retaining only the N x N matrix (1 < i,j < N). We can combine their 



technique with the reparametrization identities. Therefore we obtain a nonhnear 
renormahzation group equation in the case of a single couphng constant 



('V^ + 2)F(iV.«) = G(9,|)+0(l), (38) 

f ) - -f f 4 - - - - iT - (- w - D ■ 

(39) 

where the saddle point a of a = ^n+i n+i is determined using the expectation value 
of the resolvent 

a + ga^-w(^-^^^^^ =0. (40) 

By using the same power series expansion method, we have computed the fixed point 
and the critical exponents from these equations too. We get the same result as with 
the eigenvalue method. Therefore we believe their method is equivalent to ours in 
the eigenvalue representation. 

To obtain the expression (^), we have evaluated the integration over ^n+i n+i = 
a by the saddle point method. Let us note that it is essential to perform the 
integration over a to reproduce the exact results. It is not correct to ignore the 
degree of freedom corresponding to a by the large- assumption. 

The disadvantage of the coset technique compared to ours is that one has to 
introduce auxiliary fields when interactions like ($Ar+ij $jjv+i)^ are present as in 
the case of quartic or higher degree potentials. On the other hand, our eigenvalue 
representation allows all the potentials to be treated on the same footing. Never- 
theless, the coset technique does seem to have wider applicability, especially if we 
try to consider the cases unsolved so far. Therefore we would like to consider this 
method as a possible convenient technique for examining matrix models such as 
those corresponding to c > 1. 

From our study we have now learned that the linearized approximation neglecting 
the nonlinear terms in our renormahzation group equation provides a result not too 
far from the exact result, at least in the cases which we have studied. Therefore we 
hope to employ the linearized approximation to study models admitting no exact 
solution. This is equivalent to ignoring higher order terms with respect to the 
coupling constant g in the right hand side G of the renormahzation group equations. 
We stress again that this approximation should be done only after reparametrization 
identities are taken into account. We have already computed a number of such cases. 

The nonlinearity of the renormahzation group equations also makes it difficult 
to draw the renormahzation group flow or to characterize the m = 2 critical line as 
a trajectory of the flow. However, once we make a linearized approximation, it is 
easy to know the flow vector field in the coupling constant space. Therefore we can 
use the linearized renormahzation group equation to visualize the renormahzation 
group fiow approximately or qualitatively. This approximation should be good near 
the origin. We are currently trying to find a scheme for a more precise study of the 
renormahzation group fiow. 



It is straightforward to write down the renormahzation group equation and the 
saddle point equation for any type of matrix model It would be nice to reduce 
the general multi-matrix reparametrization identities (loop equations) to the more 
tractable form of a single matrix case. We have already observed that our method 
can be extended in a straightforward way to deal with the two-matrix model, which is 
known to describe any minimal matter coupled to two-dimensional gravity. We can 
also find operator content at any fixed points in our approach, as in eqs. (|30D - ( p2|) 
of the one-matrix case. We hope to report the results in subsequent publications. 

Note added: After we finish writing up this paper, we noticed that a preprint 



13| has appeared. In the paper, he analyzed the one-matrix model with the cubic 
potential using the coset technique. He used the first several reparametrization 
identities, but has not been able to take account of contributions from all higher 
induced interactions which are of the same order. He has not included the a-integral. 
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authors (S.H.) thanks J. Zinn- Justin, I. K. Kostov, E. Brezin, and P. Di Vecchia, for 
discussion and hospitality during his visit to Saclay, Ecole Normale, and Nordita. 
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stage of this work. This work is supported in part by Grant-in-Aid for Scientific 
Research (S.H.) and (No.05640334)(N.S.), and Grant-in-Aid for Scientific Research 
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